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We carry out a systematic study of the static and dynamic polarizabilities of the potassium atom 
using a first-principles high-precision relativistic all-order method in which all single, double, and 
partial triple excitations of the Dirac-Fock wave functions are included to all orders of perturbation 
theory. Recommended values are provided for a large number of electric-dipole matrix elements. 
Static polarizabilities of the 4s, <ipj, 5s, 5pj, and 3dj states are compared with other theory and 
experiment where available. We use the results of the polarizability calculations to identify magic 
wavelengths for the 4s — np transitions for n — 4, 5, i.e. those wavelengths for which the two levels 
have the same ac Stark shifts. These facilitate state-insensitive optical cooling and trapping. The 
magic wavelengths for the 4s — 5p transitions are of particular interest for attaining a quantum gas 
of potassium at high phase-space density. We find 20 such wavelengths in the technically interest 
region of 1050 — 1130 nm. Uncertainties of all recommended values are estimated. 

PACS numbers: 31.15.ac, 37.10.Dc, 31.15.ap, 31.15.bw 



I. INTRODUCTION 

Due to their applications in ultra-precise atomic clocks, 
degenerate quantum gases and quantum information, the 
magic wavelengths of atoms have become a subject of 
great interest in both experiments and theory 0- 

Il5| . The energy levels of atoms trapped in a light field 
are shifted by an amount that is proportional to their 
frequency-dependent polarizability, so the difference in 
the energies of any two levels depends upon the trapping 
field. This difference is often called the "ac Stark shift" . 

The idea of a "magic" wavelength, A mag i C , at which 
there is no relative shift of a give n pair of energy lev- 
els, was first proposed in Refs. [lja LLZl in the context of 
optical atomic clocks. An atom confined in a trap con- 
structed of light with a magic wavelength for the clock 
transition will, to lowest order, have the same transition 
energy as it does in free space. 

This simple idea has a number of other applications. 
A problem arises in cooling and trapping schemes, where 
the ac Stark shift of the cooling or trapping transition 
may lead to heating. Recent experiments in 6 Li [l8j and 
40 K Q degenerate quantum gases in optical traps demon- 
strated temperature reductions by a factor of about five 
and phase-space density increases by at least a factor of 
ten by laser cooling using ultraviolet (UV) transitions 
(2s — 3p and 4s — 5p, respectively) compared to con- 
ventional cooling with the visible or infrared D\ and D 2 
transitions. However, the ac Stark shifts due to trap light 
must be nearly the same for both levels in the transition 
to allow for efficient and uniform cooling . This is ac- 
complished by building the optical trap using light with 
the magic wavelength for the corresponding UV transi- 
tions. The use of the magic wavelengths is also advanta- 



geous for trapping and controlling atoms in high-Q cavi- 
ties in the strong coupling regime, so as to minimize de- 
coherence in quantum computation and communication 
protocols [lj| , and in the implementation of the Rydberg 
gate for quantum computing with neutral atoms [2(J [21j . 

Variations on the magic wavelength idea include the 
use of multiple light fields to attain ac Stark shift targets 
22] or to maximize differential response between different 
atomic species - for example, the "tune-out" wavelengths 
that trap one species but not another 0, [24| . Design and 
evaluation of all these applications requires accurate data 
on atomic frequency-dependent polarizabilities. One goal 
of our present work is to provide a list of all magic wave- 
lengths for potassium UV 4s — 5pj transitions in regions 
that are convenient for laser cooling of ultracold gases 
to high phase-space densities. For example, in 2011, 
low-temperature high-density magneto-optical trapping 
of potassium using the open 4si/2 — > 5p3/2_transition at 
405 nm was performed by McKay et al. f2|. Fermionic 
40 K was captured using a magneto-optical trap (MOT) 
on the closed 4s — > 4p 3 / 2 transition at 767 nm and then 
transferred, with high efficiency, to a MOT on the open 
4s — > 5p 3 /2 transition at 405 nm. Because the 5p 3 / 2 state 
has a smaller linewidth than the 4p 3 / 2 state, the Doppler 
limit is reduced from 145/iK to 24/iK, and temperatures 
as low as 63(6)/iK were observed. 

In this paper we provide a list of magic wavelengths 
for the 4s — 4p and 4s — 5p transitions, calculate dc 
and ac polarizabilities for several low-lying states, and 
provide recommended values for a number of relevant 
electric-dipole transitions which are of interest to appli- 
cations such as those described above. Where possible, 
we compare our results with available experimental (25| 
and high-precision theoretical values [26j . 
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Some of the calculations reported here required evalua- 
tion of the electric-dipole matrix elements for very highly 
excited states, such as 14s. These states are needed since 
the ac polarizabilities for the magic wavelengths of par- 
ticular experimental interest (around 1050 nm) are dom- 
inated by the 5p — nl transitions with n = 12 — 14. Such 
states were previously beyond the capabilities of the all- 
order method used here due to the large spatial extent 
of the orbitals. In this work, we resolved the numerical 
problems associated with such calculations and success- 
fully demonstrated the stability of our new approach. 

We begin with a brief review of recent research on the 
applications of magic wavelength concepts in Section [II] 
The calculation of electric-dipole matrix elements, static 
and dynamic polarizabilities as well as their uncertainties 
is discussed in Section Mil The magic wavelengths are 
discussed in Section HVl 



II. REVIEW OF MAGIC WAVELENGTH 
STUDIES 



where three designs for blue-detuned dipole traps were 
presented. 

Magic wavelengths for the alkali-metal atoms from Na 
to Cs, for which the ns and npi/ 2 or np 3 / 2 atomic levels 
have the same ac Stark shifts, were evaluated by Arora 
et al. [3]. The case of circular polarization was con- 
sidered in [3(J McKeever et al. [l9| demonstrated 
state-insensitive trapping of Cs atoms at 935 nm while 
maintaining a strong coupling for the 6si/2 — Qp?,/ 2 tran- 
sition. A bichromatic scheme for state-insensitive optical 
trapping of Rb atom was explored in Ref. [22] . In the case 
of Rb, the magic wavelengths associated with monochro- 
matic trapping were sparse and relatively inconvenient. 
The bichromatic approach yielded a number of promising 
magic wavelength pairs. The precise magic wavelengths 
for Li 2s — 2pj and 2s — 3pj transitions in convenient 
wavelength regions were recently calculated in [3^ |. The 
results were presented for both 6 Li and 7 Li to illustrate 
the possibilities for differential light shifts between the 
two isotopes. 



Up to the present, most work on magic wavelengths 
has been done on group I and group II atoms, which 
are the species most easily cooled and trapped by optical 
methods. We summarize these in turn. The following 
examples are representative of significant applications of 
the magic wavelength concept, however, these examples 
are not intended to constitute an exhaustive review. 



A. Group I 

The cancellation of the differential ac Stark shift of 
the microwave hyperfine clock transition in trapped 87 Rb 
atoms was demonstrated in [3J. The technique had im- 
plications for experiments involving the precise control 
of optically trapped neutral atoms, but the cancellation 
comes at the expense of a small magnetic-field sensitivity. 
"Doubly magic" conditions in magic- wavelength trapping 
of ultracold alkali- metal atoms were investigated by Dere- 
vianko [5]. This work demonstrated that the microwave 
transitions in alkali-metal atoms may be indeed made im- 
pervious to both trapping laser intensity and fluctuations 
of magnetic fields. 

The issue of the mismatch of the polarizabilities of the 
ground and excited states has also arisen in the Ryd- 
berg gate approach to quantum information processing 
[27Ll28j] , in which the qubit is based on two ground hyper- 
fine states of neutral atoms confined in an optical lattice. 
An atom in a Rydberg state will, in general, move in a dif- 
ferent optical lattice potential than that experienced by 
the ground state. Therefore, the vibrational state of the 
atom in the lattice may change after the gate operation is 
completed, leading to decoherence due to motional heat- 
ing. This problem may be resolved by the use of magic 
wavelengths 0, Hl|. Use of magic wavelengths in op- 
tical traps for Rydberg atoms was also discussed in Q, 



B. Group II 

The magic wavelengths for the Sr 5s 2 1 <S'o — 5s5p 3 Pg 
and 5s 2 1 S — 5s5p 3 P° transitions have been measured 
in [Hill. The Yb clock 6s 2 ^0 - 6s6pfP ° magic wave- 
length was predicted to be 752 nm in [35j and measured 
to be 759.355 nm in Ref. (36|. The magic wavelength for 
the ultraviolet 6s 2 1 S'o O 6s6p 3 Pq clock transition in 
Hg was recently reported by Yi et al. [l[ . The Stark-free 
(magic) wavelength was found to be 362.53(0.21) nm, 
in excellent a gree ment with the theoretical prediction 
360 nm from [33], calculated using a method that com- 
bines configuration interaction and many-body pertur- 
bation theory. The magic wavelengths of other group II 
and group lib atoms have been predicted in [38]. The 
magic wavelengths are very sensitive to the values of the 
ac polarizabilities and allow for precise tests of the theory 
E3] • Moreover, the magic wavelengths can be used to 
determine the values of important electric-dipole matrix 
elements which are difficult to obtain by direct experi- 
mental techniques. For example, the 5s5p 3 Pq — 5s6s 3 Si 
matrix element in Sr was recently determined using the 
experimental value of the Sr 5s 2 1 Sn — 5s5p 3 P£ magic 
wavelength with 0.5% precision [40(. Dammalapati et 
al. @ investigated light shifts of heavy alkaline earth 
elements barium (Ba) and radium (Ra) , which are of in- 
terest for development of optical lattice clocks and for 
permanent electric dipole moment searches. The wave- 
length dependence of light shifts of the ns 2 1 Sq ground 
state, the nsnp 3 P 1 ° and ns(n — l)d 1 D 2 excited states in 
Ba (n = 6) and the ns 2 1 S'o ground state, the nsnp 3 P 1 ° 
and ns(n — l)d 3 D 2 excited states in Ra (n = 7) were 
calculated. Several magic wavelengths in the visible and 
infrared regions accessible with commercial lasers for op- 
tical dipole trapping of Ba and Ra were identified [9(. 
Magic wavelengths of an optical clock transition of bar- 
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TABLE I: Absolute values of the reduced electric-dipole matrix elements in K and their uncertainties in a.u.. The present 
all-order values are given unless noted otherwise. ' a 'Expt. [2^|, ^'determined from Stark shift data in [lj. The uncertainties 
are estimated where possible (see text). 
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ium were presented in 10]. Dipolc polarizabilities of 
ns 2 1 Sq and nsnp 3 Pq states and relevant magic wave- 
lengths of Sr, Yb, Zn, Cd, and Hg atoms were studied by 
a semiempirical approach in Refs. [13. [T^. 

The magic wavelength conditions that can make opti- 
cal lattice clocks insensitive to atomic motion were pre- 
sented by Katori et al. 6]. This work demonstrated 
that the spatial mismatch of the interactions in the clock 
transition can be treated as a spatially constant offset 5v 
for specific lattice geometries. Numerical estimates were 
made for Sr @. 

Theoretical study of the dynamic scalar polarizabili- 
ties of the ground and selected long-lived excited states 
of dysprosium was recently carried out by Dzuba et al. 
• A set of the magic wavelengths of the unpolarized lat- 
tice laser field for each pair of states, which includes the 
ground state and one of these excited states was given. 
The authors presented an analytical formula that allows 



for the determination of approximate values of the magic 
wavelengths without calculating the dynamic polarizabil- 
ities of the excited states near resonances I7J . 



III. MATRIX ELEMENTS AND 
POLARIZABILITIES 

The magic wavelengths for a specific transition are 
located by calculating the frequency-dependent polar- 
izabilities of the lower and upper states and finding 
their crossing points. The all-order approach to the 
calculation of atom ic p olarizabilities was discussed in 
Refs. 0, H [H, |H, S3 S3, and we provide only a brief 
summary of the methods here. Unless stated otherwise, 
all specific data refers to the K atom, and we use the 
conventional system of atomic units, a.u., in which e, m e , 
47reo and the reduced Planck constant Ti have the numer- 
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TABLE II: Values of scalar (ao) and tensor (02) polarizabil- 
ities in K. The present results are compared with theoretical 
and experimental values. Ref. [26j did not include uncertainty 
estimates. All values are in atomic units. 
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ical value 1. Polarizability in a.u. has the dimension of 
volume, and its numerical values presented here are ex- 
pressed in units of ag, where ao « 0.052918 nm is the 
Bohr radius. The atomic units for a can be converted 
to SI units via a/h [Hz/(V/m) 2 ]=2.48832xl0- 8 a [a.u.], 
where the conversion coefficient is 47reoao//i and the 
Planck constant h is factored out. 

The frequency-dependent scalar polarizability, a(uj), of 
an alkali-metal atom in the state v may be separated into 
a contribution from the ionic core, a corc , a core polariz- 
ability modification due to the valence electron, a vc , and 
a contribution from the valence electron, a v {uj). The va- 
lence contribution to frequency-dependent scalar ao and 
tensor a?2 polarizabilities is evaluated as the sum over 
intermediate k states allowed by the electric-dipole tran- 
sition rules [ifl 

V( s = 2 {k \\ d \\ v f( Ek -E v ) 
aa[UJ > 3(2^ + 1)^ {E k -E v y-^ ' 

osm = -icY,(~^ +3k+1 [ J ; I \ } 

(E k - E v y - a; 2 ' { ' 

where C is given by 

c _( 5j„(2j„-l) \ 1/2 

v 60„ + l)(2^ + l)(2j t , + 3)y' 

and (A: ||d|| u) are the reduced electric-dipole matrix ele- 
ments. In these equations, uj is assumed to be at least 
several linewidths off resonance with the corresponding 
transitions. Linear polarization is assumed in all calcula- 
tion. The ionic core polarizability and a vc term depend 
weakly on uj for the frequencies treated here and are ap- 
proximated by their dc values calculated in the random- 



phase approximation (RPA) (45[. We find the contribu- 
tion from the K + ionic core to be a C orc = 5.457 clq. A 
counter term a vc compensating for excitation from the 
core to the valence shell which violates the Pauli princi- 
ple is very small. For example, it is a vc = —0.00015 a.u. 
for the 4pj states of K. 

We use the linearized version of the coupled clus- 
ter approach (also referred to as the all-order method), 
which sums infinite sets of many-body perturbation the- 
ory terms, for all significant terms in the equations above. 
The 4s— np, Ap—nl, 5s— nl, 5p— nl, and Zd—nl transitions 
with n < 26 are calculated using this approach 0, |46| . 

As we noted in the Introduction, the present calcu- 
lation required evaluation of the electric-dipole matrix 
elements for highly excited states, since the frequency- 
dependent polarizabilities for the 4s — 5p magic wave- 
lengths of particular experimental interest are dominated 
by the 5p — nl transitions with n = 12 — 14. The diffi- 
culty with the applications of the all-order method for 
these states results from the use a complete set of Dirac- 
Fock (DF) wave functions on a nonlinear grid gener- 
ated using B-splines constrained to a spherical cavity. 
A large cavity with radius of R = 220 ao is needed to 
accommodate all valence orbitals with ns = 4s — 10s, 
np = Ap — IQp, and nd = 3d — 9d. A cavity radius of 
400 ao was chosen to accommodate additional valence 
orbitals with ns = lis — 14s, np = lip — I3p, and 
nd = lOd — I2d. Our basis set consists of 70 splines of or- 
der 11 for each value of the relativistic angular quantum 
number k for R — 220 ao and 100 splines of order 13 for 
R = 400 a Q . We have conducted test comparisons of the 
basis set energies with the actual DF values to demon- 
strate the numerical stability of this calculation. We can 
use available experimental energies for the ns = 4s — lis, 
np = Ap — lOp, and nd = 3d — I2d states from [13] and 
theoretical all-order energies for other states with n < 26. 
The remaining small contributions with n > 26 are cal- 
culated in the DF approximation. For example, the con- 
tributions from states with n > 26 give only 0.075 a.u. 
to the polarizability of the Ap 1 / 2 state. We note that 
states with n > 19 in our basis have positive energies 
and provide a discrete representation of the continuum. 

The evaluation of the uncertainty of the matrix ele- 
ments in this approach was described in detail in [32], [44| ■ 
Four all-order calculations were carried out, including 
two ab initio all-order calculations with and without 
the inclusion of the partial triple excitations and two 
other calculations that incorporated semiempirical esti- 
mates of high-order correlation corrections starting from 
both ab initio runs. The spread of these four values 
for each transition defines the estimated uncertainty in 
the final results when considered justified based on the 
dominant correlation contributions to the El matrix el- 
ements [H], HH . We note that this procedure does not 
work in the small number of cases where we can not 
estimate uncertainty in the dominant contributions us- 
ing the procedure described above. No uncertainties are 
listed for such matrix elements, however, their contribu- 



5 

9s 7d 5/2 8s 6d 5/2 7s 5d 5/2 6s 4d 5/2 4s 




500 550 600 650 700 750 800 

Wavelength X in vacuum (nm) 



FIG. 1: (Color online) The frequency-dependent polarizabilities of the K 4s and 4p 3 / 2 states. The magic wavelengths are 
marked with circles. The approximate positions of the 4p 3/ / 2 — nl resonances are indicated by vertical lines with small arrows 
on top of the graph, together with the corresponding nl. 



tions are small, leading to negligible effects on the final 
uncertainties of the polarizabilities. The absolute val- 
ues of the reduced electric-dipole matrix elements used 
in our subsequent calculations and their uncertainties are 
listed in a.u. in Table HI We list only the most important 
subset of the several hundred matrix elements that were 
calculated in this work. 

Our results for scalar and tensor polarizabilities of the 
Apj excited states of potassium are compared with calcu- 
lations of [48l l49l j and with experimental measurements 
reported by Marrus and Yellin [42[ in Table|n] The Bates- 
Damgaard method was used by Schmieder et al. [48| and 
the time-dependent gauge-invariant variational method 
was used by Marrus and Yellin [42f . The uncertainty 
in the experimental measurement [42j of the scalar po- 
larizability is too large to reflect on the accuracy of the 
present calculations. Extensive comparison of the the- 
oretical and experimental static polarizabilities for the 
alkali- metal atoms was recently given in the review [45j . 



IV. MAGIC WAVELENGTHS 



14s 12d 




1051.5 1052.0 1052.5 1053.0 1053.5 1054.0 1054.5 1055.0 
Wavelength \ in vacuum (nm) 
FIG. 2: (Color online) The frequency-dependent polarizabili- 
ties of the K 4s and 5p 3 / 2 states. The magic wavelengths are 
marked with circles and arrows. The approximate positions 
of the 5p 3 / 2 — 14s and 5p 3 / 2 — 12d resonances are indicated 
by vertical lines with small arrows on top of the graph. 



We define the magic wavelength A mag i c as the wave- 
length for which the ac polarizabilities of two states in- 
volved in the atomic transition are the same, leading to 
a vanishing ac Stark shift of that transition. For the 
ns — np transitions, a magic wavelength is represented 
by the point at which two curves, a ns (X) and a np (\), 



intersect as a function of the wavelength A. The total 
polarizability for the np 3 / 2 states is given by a = ao — a i 
for rrij = ±1/2 and a — ao + a 2 for the rrij = ±3/2 case. 
Therefore, the total polarizability of the np 3 / 2 state de- 
pends upon its rrij quantum number and the magic wave- 
lengths need to be determined separately for the cases 
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FIG. 3: (Color online) The frequency-dependent polarizabilities of the K As and 5p 1 /2 states. The magic wavelengths are 
marked with circles. The approximate positions of the 5pi/2 — nl resonances are indicated by vertical lines with small arrows 
on top of the graph, together with the corresponding nl. 



14s 12d 3/2 13s lld 3 / 2 12s 10d 3/2 lis 9d 3/2 




1040 1060 1080 1100 1120 1140 



Wavelength X in vacuum (nm) 



FIG. 4: (Color online) The frequency-dependent polarizabilities of the K 4s and 5p3/2 states. The magic wavelengths are 
marked with circles. The approximate positions of the 5p%/2 — nl resonances are indicated by vertical lines with small arrows 
on top of the graph, together with the corresponding nl. 
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TABLE III: Magic wavelengths for the 4s — npj transitions 
in K. The 500 - 1227 nm and 1050 - 1130 wavelength ranges 
were considered for the 4s — 4pj and 4s — 5pj transitions, 
respectively. The corresponding polarizabilities are given in 
a.u. The resonance near the magic wavelengths are listed in 
the first column. 

Resonance K A ma gic a 



with nij — ±1/2 and rrij = ±3/2 for the ns — np^n tran- 
sitions, owing to the presence of the tensor contribution 
to the total polarizability of the np 3 / 2 state. The uncer- 
tainties in the values of magic wavelengths are found as 
the maximum differences between the central value and 
the crossings of the a ns ± Sa ns and a np ± 5a np curves, 
where the 8a are the uncertainties in the correspond- 
ing ns and np polarizability values. All calculations are 
carried out for linear polarization. Several magic wave- 
lengths were calculated for the 4s — 4p!/ 2 and 4s — 4p 3 / 2 
transitions in K in Ref. 14j] using the all-order approach. 
Only the magic wavelengths with A > 600 nm were listed. 
In this work, we present several other magic wavelengths 
for these D\ 1 D 2 transitions above 500 nm. 

The frequency-dependent polarizabilities of the 4s and 
4^3/2 states for A = 500 — 800 nm are plotted in Fig. Q] 
The magic wavelengths are marked with circles. The 
approximate positions of the 4p 3 / 2 — nl resonances are 
indicated by vertical lines with small arrows on top of 
the graph, together with the corresponding nl. For ex- 
ample, the arrow labelled 7s indicates the position of 
the 4p 3 / 2 — 7s resonance. The corresponding magic 
wavelengths are listed in Table IIIII We note that the 
4j>3/2 — 5s resonance wavelength is outside of the plot 
region at A = 1253 nm). While there are 8 magic wave- 
lengths for the 4s — 4p 3 / 2 \nij\ — 1/2 transition in the 
wavelength region shown on the plot, there are only 4 
magic wavelengths for the 4s — 4p 3 / 2 \mj\ = 3/2 tran- 
sition since there are no corresponding crossings near 
the 4p 3 / 2 — ns resonances as in the case of \nrij\ = 1/2. 
The 769 nm magic wavelength for the \m,j\ — 1/2 is 
not shown on the plot since the corresponding polariz- 
ability (-27400 a.u.) is outside of the plot y-axis range. 
There is only one magic wavelength above 800 nm for 
the 4s — 4p 3 / 2 \m,j\ = 1/2 transition due to 4p 3 / 2 — 5s 
resonance and none for the \nij\ = 3/2 case. The magic 
wavelengths for the 4s — ^Pi/2 transition are very close 
to those for 4s — 4p 3 / 2 \nij\ — 1/2. They are also given 
in TableUm 

The magic wavelengths for the UV 4s — 5pj transi- 
tions are completely different than those for the D±, £> 2 
lines owing to completely different set of resonances. The 
K case is also significantly different from that of Li [32j 
due to differences in the resonant transition wavelengths. 
We list the magic wavelengths for the 4s — 5pi/2 and 
4s— 5p3/2 transitions in the range of 1050-1130 nm, which 
is of particular experimental interest in Table [TTTl We find 
20 magic wavelengths in the technically interest region of 
1050 — 1130 nm accessible by a number of widely used 
lasers. The magic wavelengths for the 4s — 5p 3 / 2 transi- 
tion near 1053 nm wavelength are illustrated in Fig.[2j As 
in the case of the 4s — 4p 3 / 2 transition, there is no magic 
wavelength for the \m,j \ — 1/2 case near the ns resonance. 
All magic wavelengths for the 4s — 5pi/ 2 and 4s — 5p 3 / 2 
transitions in the range of 1050-1140 nm are illustrated 
in Figs. [3] and Figs.0J The same designations are used as 
in the previous graphs. Comparing these figures with the 
similar plots for Li (see Figs. 3 and 4 of Ref. [12]) shows 
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that K magic wavelengths near 1050-1130 nm originate 
from crossings near much higher resonances (n = 9 — 14 
vs. n = 6 — 7 for Li) making the calculation for K more 
complicated due to the very large cavity size required to 
accommodate such highly excited orbitals. 

V. CONCLUSION 

We have calculated the ground 4s, 4p, and 5p state 
ac polarizabilities in K using the relativistic linearized 
coupled-cluster method and evaluated the uncertainties 
of these values. We have used our calculations to iden- 



tify the magic wavelengths for the 4s — 4p and 4s — 5p 
transitions. The magic wavelengths for the ultraviolet 
resonance lines is of particular interest for laser cooling 
of ultracold gases with high phase-space densities. 
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